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Abstract 

The relativistic heavy ion program developed at RHIC and now at LHC motivated a deeper study 
of the properties of the quark gluon plasma (QGP) and, in particular, the study of perturbations 
in this kind of plasma. We are interested on the time evolution of perturbations in the baryon 
and energy densities. If a localized pulse in baryon density could propagate throughout the QGP 
for long distances preserving its shape and without loosing localization, this could have interesting 
consequences for relativistic heavy ion physics and for astrophysics. A mathematical way to proove 
that this can happen is to derive (under certain conditions) from the hydrodynamical equations 
of the QGP a Korteveg-de Vries (KdV) equation. The solution of this equation describes the 
propagation of a KdV soliton. The derivation of the KdV equation depends crucially on the 
equation of state (EOS) of the QGP. The use of the simple MIT bag model EOS does not lead to 
KdV solitons. Recently we have developed an EOS for the QGP which includes both perturbative 
and non-perturbative corrections to the MIT one and is still simple enough to allow for analitycal 
manipulations. With this EOS we were able to derive a KdV equation for the cold QGP. 



I. INTRODUCTION 



Korteweg - de Vries solitons are very interesting non-linear waves, which may exist in 
many types of fluids from ordinary water to astrophysical plasmas In the last years we 
have started to produce a new kind of fluid in laboratory: the quark gluon plasma (QGP). 
This is a state where quarks and gluons, usually confined in the interior of baryons (such 
as the proton) and mesons, are free to travel longer distances. With the beginning of the 
LHC era, we have means to study larger and longer living samples of QGP and even the 
propagation of perturbations in this new medium. In this context a natural question is: can 
we have KdV solitons in the QCD plasma? In this work we give an answer to this question. 

Before the QGP there were other fluids made of strongly interacting hadronic matter and 
the existence of KdV solitons in these fluids was already investigated. The first works on the 
subject were published in [2], where the authors considered the propagation of baryon density 
pulses in proton-nucleus collisions at intermediate energies. In this scenario the incoming 
proton would be absorbed by the nuclear fluid generating a KdV soliton, which, traversing 
the whole nucleus without distortion, would escape from the target as a proton and would 
simulate an unexpected transparency. In 2] the existence of the KdV soliton relied solely 



on the equation of state (EOS), which had no deep justification. In [3( we have reconsidered 
the problem, introducing an equation of state derived from relativistic mean field models 
of nuclear matter. We concluded that the homogeneous meson field approximation was 
too strong and would exclude the existence of KdV solitons. We could also trace back 
the derivative terms in the energy density to derivative couplings between the nucleon and 
;he vector meson. In j^J we extended our analysis to relativistic hydrodynamics and in 
5] to spherical and cylindrical geometries. In p we considered hadronic matter at finite 
temperature and studied the effects of temperature on the KdV soliton. In [7| we started the 
study of perturbations in the QGP at zero and finite temperature. The conclusion found in 
that work was that the existence of KdV solitons in a QGP depends on details of the EOS 
and with a simple MIT bag model EOS there is no KdV soliton! A further study of the 
equation of state, carried out in sj], showed that if non-perturbative effects are included in 
the EOS through gluon condensates, then new terms appear in the expression of the energy 
density and pressure and in the present work we show how these new terms lead to a KdV 
equation, after the proper treatment of the hydrodynamical equations. 
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In the next section we briefly review the equations of one- dimensional relativistic fluid 
dynamics. In section III we introduce the equation of state, in section IV we derive the KdV 
equation and in section V we present a numerical analysis of the obtained equation. 



II. RELATIVISTIC FLUID DYNAMICS 



Relativistic hydrodynamics is well presented in the textbooks 
version of the Euler equation 7], |9_|, [l0| is given by: 



dv 
dt 



+ (v-V)v 



(,w Vp+l7 



dp 
dt 



The relativistic 



(1) 



where v, e, p and 7 are the velocity, energy density, pressure and the Lorentz factor respec- 
tively. We employ the natural units c = 1 and h = 1. Space and time coordinates will be 
in fm (lfm = 10 _15 m). The relativistic version of the continuity equation for the baryon 



density is 



9 : 



duj B v = 

Since js v = u l ' Pb the above equation can be rewritten as Q: 

dps 2 ^ f dv _ -*\ -l . . 
+7 v p B I + v ■ Vv I + V - [pb v) 







(2) 



(3) 



where ps is the baryon density. In the one dimensional Cartesian relativistic fluid dynamics 
the velocity field is written as v = v(x,t)x where x is the unit vector in the x direction. 
Equations ([1]) and ([3]) can be rewritten in the simple form: 



dv dv (f 2 — I) ( dp dp\ 
(e + p) \dx dt J 



dt ~^ V 'dx 



and 



vpB 



f dv dv\ 

Kbl + v d-x) 



, n 2^( d PB . dv dp B 
+ ^- V \-dT + pB d-x +V ^x- 



(4) 



(5) 



III. THE QGP EQUATION OF STATE 



In what follows we present 



works on the subject see 



the mean field treatment of QCD developed in [8| (for previous 
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12( | ) and go beyond the homogeneous field approximation, 



including the terms with gradients. 
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The Lagrangian density of QCD is given by: 

1 Nf - r 

Cqcd = —F; V F^ + £ # [iT M (%^ - igT?pi) - 5 %3 m q ] ^ (6) 

4 9=1 

where 

= d^G au - d v G a » + gf abc G b »G cv (7) 

The summation on g runs over all quark flavors, m q is the mass of the quark of flavor q, 
i and j are the color indices of the quarks, T a are the SU(3) generators and f abc are the 
SU(3) antisymmetric structure constants. For simplicity we will consider massless quarks, 
i.e. m„ = 0. Moreover, we will drop the summation and consider only one flavor. At the 

nn 

end of our calculation the number of flavors will be recovered. Following [111 . 1121 ]. we shall 
write the gluon field as: 

where A a ^ and a a ^ are the low ("soft") and high ("hard") momentum components of the 
gluon field respectively. We will assume that A ail represents the soft modes which populate 
the vacuum and the terms containing A a>1 will be replaced by their expectation values (v4 a/i ), 
(A aM Ap, etc... in the plasma. a a>1 represents the modes for which the running coupling 
constant is small. 

In a cold quark gluon plasma the density is much larger than the ordinary nuclear matter 
density. These high densities imply a very large number of sources of the gluon field. Assum- 
ing that the coupling constant is not very small, the existence of intense sources implies that 
the bosonic fields tend to have large occupation numbers at all energy levels, and therefore 
they can be treated as classical fields. This is the famous approximation for bosonic fields 
used in relativistic mean field models of nuclear matter [k| • It has been applied to QCD in 
the past and amounts to assume that the "hard" gluon field, a^, is simply a function of the 
coordinates: 

a%(x,t) = 8 lja a%(x,t) (9) 
with d v a a 7^ 0. This space and time dependence goes beyond the standard mean field 



approximation 



13|, where a a is constant in space and time and consequently d v a a = 0. We 



keep assuming, as in |8[, that the soft gluon field A a>Ji is independent of position and time 
and thus d v A aii = . Following the same steps introduced in j^] we obtain the following 
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effective Lagrangian: 



£o = -2«o( v «o) + -^-"o"o - B QCD 



+ ^ i (i5 ljl »d, + g 1 Q T«ai)^ j 



(10) 



where mc is the dynamical mass of the hard gluon a generated by its interaction with 
the soft gluons A a ^ from the vacuum and it is related to the dimension two (A 2 ) gluon 
condensate. The constant Bqcd is related to the dimension four gluon condensate (F 2 ) (see 
|s | for details). 

The effective Lagrangian ( TTUj) is quite similar to the one obtained in jg] and the only 
difference is the first term, which is new and comes from the gradients. The equations of 
motion jsj are given by: 

dC 



— -d 9C I 8 d 
drji M <9(<V?;) v " 



(11) 



Inserting ( ITU]) into (fTTj) with r\\ = ocq(x, t) and r] 2 = if)(x,t) we find: 

- V 2 a a Q + m G 2 a% = -gp a (12) 

(tYd, + 91 °T a a^ = (13) 

where p a is the temporal component of the color vector current given by j au = ^i r fTf^)j . 
The energy-momentum tensor reads {(J: 



dC 



dC 



dC 



d(d^dprii 



(14) 



From the above expression we can obtain the energy density (e —< T 00 >) which turns out 
to be bl: 



- = ^«o(V 2 « a ) 



iUG -a a a a + B QCD -gp a a a + 3? Q 



dkk V k 2 + m 2 



(15) 



2 " " ' •" " ' 2ir 2 Jo 

where jq is the quark degeneracy factor 7q = 2(spin) x 3(flavor). The sum over all the color 
states was already performed and resulted in the pre-factor 3 in the expression above. k F is 
the Fermi momentum defined by the quark number density p: 

p = (N\4^\N) = 1 <£(N\N) = / d*k = 3^ f*' dk k 2 



A, A 



lQ :k/ 



2tt 2 



(16) 



In the above expression \N) denotes a state with N quarks. In a first approximation the 
ield «q rnay be estimated from ( TT2l) . Neglecting the derivative term V 2 «o of ( [121) we have: 



Inserting ( ITT)) in the first term of (fT2l) and then solving it for «q we find: 

771(3 TJ-G 

We can write the color charge density p a in terms of the quark number density p through: 

p a p a = 3p 2 (19) 

Analogously we have 

p a VV = 3pV 2 p, p a V 4 p a = 3pV 4 p (20) 

Inserting (ITS)) . (IT9)) and (120)) into (ITS)) , performing the momentum integral and using the 
baryon density, which is p B = \p, we arrive at the final expression for the energy density in 
one spatial dimension: 

_ ( Kg 2 \ 2 / 27g 2 \ d 2 p B / 27g 2 \ d 4 p B f 21 g 2 \ d 2 p B d 4 p B 
E \2m G 2 J pB \2m G A ) PB dx 2 \2m G & ) pB dx 4 \2m G 8 J dx 2 dx A 

The pressure is given by p = | < Ta >. Repeating the same steps mentioned before we 
arrive at: 



V 



( 27 'g 2 \ 2 + rl8g 2 \ d 2 p B _ / 9g* \ &p B _ / 9g 2 \dp B dp B 
\2m G 2 J \m G 4 J dx 2 \m G 6 J dx 4 \2m G 4: J dx dx 



| / 9g 2 \ d 2 p B d 2 p B _ / 9g 2 \ d 2 p B d 4 p B _ / 9g 2 \d 3 p B d 3 p B _ / 9g*_\ dps (Pps 
\2m G 6 J dx 2 dx 2 \m G 8 J dx 2 dx A \2m G 8 J dx 3 dx 3 \m G 6 J dx dx 3 

+ m 

where now k,F defined by ( tTBl) is given by p B = kp^/n 2 . 

IV. THE KDV EQUATION 

We now combine the equations @ and 05j to obtain the KdV equation which governs 
the space-time evolution of the perturbation in the baryon density. We first write (HI) and 
(|5j in terms of the dimensionless variables: 

P= — , v = — (23) 
Po c s 
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where po is an equilibrium (or reference) density, upon which perturbations may be gene- 
rated, and c s is the speed of sound. Next, we introduce the £ and r "stretched" coordinates 



12 {x — C s t) 



I Cat 



r = a 



(24) 



3/2 ; 

R ' R 

where a is a small expansion parameter, R is a typical size scale of the problem. After this 
change of variables we expand f l23l) as: 



p — 1 + O-pi + cr p 2 + . . . 
,2, 



(25) 

v = av\ + a^f 2 + . . . (26) 

Having rewritten (j3J) and (jSJ) in the £ — r space and having expanded them in powers of a 
up to <t 2 we organize the two equations as series in powers of a. After these steps 01]) and 
(J5J) become: 



+ 



27<? 2 p 2 



27s 2 Po 2 



+ 



27# 2 p : 



+ vr 2 / 3 po 4/3 



Y27s 2 p 2 



c s 2 + 3vr 2 / 3 p 4/3 c s 2 



+ 7T 2 / 3 p 4 / 3 



5p2 

9e 



27^ 2 p s 



m G 



c s 2 + 3tt 2 / 3 p 4 / 3 c s 2 



dpi] 
d£ j 



d^ +vl ^ 



+ 



/27s 2 po 2 



m G 



)2 Cs 2 + 4vr 2 / 3 p 4 / 3 c s 2 



Pi- 



9f i 



+ 



<9£ 

i8g 2 p 2 \a 3 p 

m G 4 R 2 J d£ 3 



V m G 2 
27ff 2 p 2 
m G 2 



Pi% + vr 2 /V /3 ^^ 



<9£ 

2/3 n 4/3^ 2 



C/ + 7T^ J p 



3 <9£ 
dpi 



(27) 



and 



9f i <9p 



<9t> 2 dp 2 dpi 



dvi 



dpi 



dvi 



I uui b-W 2 I . ^fl . ""I . "PI 2 UL >1 ( n fo«A 

respectively. In the last two equations each bracket must vanish independently and so 
{. . .} = 0. From the first term of ( 128|) we obtain pi = v\. Using this identity in the first 
term of ( )27|) we obtain an equation, which solved for c s yields: 

+ vr 2 / 3 p 4 / 3 



27g" po 

2 



Zffi) + 3vr2/3 po 4/3 



(29) 



Inserting these results into the terms proportional to a 2 , we find, after some algebra the 
KdV equation: 



dpi 
dr 



+ 



(2-c s 2 ) /27^ 2 p 2 \(2c s 2 -l) 



V m G 2 



2A 



l Cs "6 



Pi 



dpi 
3£ 



+ 



- 9fl 2 Po 2 - 
m G A R 2 A_ 



d 3 pi 
dC 3 



where: 



A 



f27g 2 p < 



V m G 2 

Returning to the x — t space we obtain: 

dpi dpi 
h c s 



+ 37T 2 / 3 P0 4/3 C S 



at 



where 



and 



a: 



{2-c s 



dx 

P = 



ac s pi 



(2c, 



dp 





\ m G 2 
d 3 pi 



+ 7T 2 / 3 P 4/3 



dx 







2A 

% 2 Po 2 c s 



dx 3 
7r 2 / 3 p 4 / 3 



(30) 



(31) 



(32) 



(33) 



(34) 



m G 4 A 

which is the KdV equation at zero temperature for the small perturbation in the baryon 
density p\ = op\. If we neglect the space and time derivatives in ( 12TT) and ( 1221) and repeat 
the derivation sketched above we arrive at: 



dpi dpi A dpi 

at ox ox 







(35) 



which is a breaking wave (BW) equation for p\. We close this section emphasizing that Eq. 
( )32|) is the main result of this work. It shows that for suitable choices of the parameters a 
and j3 we can have KdV solitons in a quark gluon plasma. 



V. NUMERICAL ANALYSIS 

The KdV equation (|32|) has an analytical soliton solution given by [lj]: 
Pi(x,t) 



3(u — c. 



sech 



[u - c s 



4/3 



(x — ut) 



(36) 



where u is an arbitrary supersonic velocity. A soliton is a localized pulse which propagates 
without change in shape and in this case it has the width A defined by: 



A 



4/3 



[u 



3Qg 2 p 2 c 



c s ) \ (u - c s )m G 4 A 



(37) 



For the numerical estimates we shall use the following values of the parameters : Bq G d = 
0.0006 GeV 4 , g = 0.35, m G = 290 MeV and p = 2fm~ 3 which, when sustituted in (EHD 



yield c s 2 = 0.5 (c s = 0.7). We choose u = 0.8. With the help of (ETJ) we find A = 1.7 fm for 
the soliton width and 0.6 for its amplitude. Even though this work is essentially qualitative 
the chosen numbers are well appropriate to study a realistic situation of a perturbation 
traversing the QGP. 

In Fig. 1 we show the numerical solution of ( 13 2 j) with initial condition pi(x,t = 0) and 
pi(x,t) given by (|36l) . We can observe the time evolution of the initial gaussian-like pulse 
as a well defined soliton, keeping its shape and form. This solution shows the behavior 
expected from the analytic solution. We see that using the correct input for the amplitude 
and width we obtain a pulse which propagates without distortion. This initial condition is 
a very special case of little practical interest. A real perturbation produced in the QGP will 
most likely have the "wrong" amplitude and "wrong" width. For arbitrary amplitudes the 
solution must be numerically calculated. 

In Fig. 2 we show again the numerical solution of (132]) with the initial condition given 
by ( l36l) multiplied by a factor 10. Now we observe that the initial pulse starts to develop 
secondary peaks, which are called "radiation" in the literature. Further time evolution will 
increase the strength of these peaks until the complete loss of localization. 

In Fig. 3 we show the numerical solution of (135]) with the initial condition given by (1361) . 
We observe the gradual formation of a "wall" followed by the dispersion of the initial pulse. 
In Fig. 4 the amplitude of the initial pulse (136]) is multiplied by a factor 10 and used as 
initial condition for ( 135]) . As expected the dispersion takes place much earlier than in Fig. 
3. 



VI. SUMMARY 



The main conclusion of this work is that it is indeed possible to have KdV solitons in 
QCD, provided that two conditions are satisified. The first condition is that the gluon 
field have a dynamical mass. In this case the equation of motion (Tl2l can be solved in 
the weak inhomogeneity approximation yielding ([TBI) . The existence of a dynamical gluon 
mass has been intensely discussed in the literature during the last years and seems to be 
well established (for details see the references given in |8|). For a massless gluon field we 
can only have a breaking wave equation, as it was found in our previous work Q. The 
second condition is the existence of second order derivative terms in the energy density and 
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FIG. 1: Numerical solution of (I32D with (|36p as initial condition calculated at different times. 
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FIG. 2: The same as Fig. 1 with an initial amplitude ten times larger. 



pressure. These terms appear naturally from the formalism, as we can see in (fT5l) . However it 
is necessary to keep these derivative terms. The use of uniform field approximations prevents 
us from finding KdV solitons. If we neglect the derivatives we arrive at the breaking wave 
equation (135|) . The practical difference between perturbations governed by the KdV and 
BW equations is that the former propagate much longer keeping its localization whereas the 
latter loose localization and may generate unstable "walls" . The numerical analysis of some 
cases confirms the anticipated qualitative expectation. The application of the formalism 
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FIG. 3: Numerical solution of (I35I) with (I36I) as initial condition. 
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FIG. 4: The same as Fig. 3 with an initial amplitude ten times larger. 



developed in this work to problems in the theory of compact stars is in progress. 
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